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On-Line Data Reduction 
for Nuclear Analyzers 

Sniffing out minute amounts of radioactivity in our environment is just one of the talents 
of the versatile multichannel analyzer. Here are four multichannel analyzer systems — two 
of them brand new — that include on-line computing devices capable of anything from 
straight number crunching to completely automating the analysis. 

By Jonathan R. Cross, James A. Doub 
and John M. Stedman 



Some of the first 'computers' were really pulse- 
height analyzers. These electronic instruments were de- 
veloped by nuclear scientists to help identify and classify 
various radioactive materials. While this is still their prin- 
cipal function, they are also used for a variety of other 
purposes in medicine, electronics, environmental science, 
and many other fields. 

The contemporary pulse-height analyzer, or multi- 
channel analyzer, is much faster, has much better reso- 
lution, and can do much more than those "rooms full of 
tubes' of twenty years ago. It collects and analyzes the 
same kind of data, but it's better data, and there's more 
of it. Reducing this data — extracting all the information 
it contains — calls for complex data manipulation on a 
scale vastly greater than that of the pencil, paper, and 
slide rule methods of the early days. 

Hewlett-Packard is a manufacturer of computers, cal- 
culators, and data communications devices as well as 
multichannel analyzers, and has developed a variety of 
systems in which analyzers are supported by on-line 
data reduction facilities. Depending upon his needs, a 
user can choose a relatively inexpensive system that does 
data reduction only, or either of two more powerful sys- 
tems that offer program control of the analyzer along 
with on-line data reduction, or a super-system which 
combines the abilities of one or more analyzers with a 
digital computer. The systems arc automatic and easy to 
use. whether the user is a technician who must repeat an 
experiment many times for many different samples or a 
researcher who requires 'hands on' control of the entire 
data accumulation and reduction process. 

Four Types of Systems 

The four types of systems and their capabilities are 
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compared in the table on page 3. At the heart of all the 
systems is the HP 5401 B Multichannel Analyzer (or 
components thereof, in the case of the largest system). 
The 540 IB is a high-speed, high-resolution instrument 
(200 MHz clock rate analog-to-digital converter. SI 92 
channels) with very low integral nonlinearity (less than 
±0.05%). Readers who aren't familiar with what mul- 
tichannel analyzers do and how they're used will find 
more information about them on pages 4 and 6. 
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1 Cover: This is one of the 
reactors at General Electric 
Company's Vallecitos Nu- 
clear Center. One of the 
Hewlett-Packard multichan- 
nel analyzer systems de- 
scribed in this issue is 
used at Vallecitos for nu- 
clear research and effluent 
monitoring. We thank GE for allowing us to pho- 
tograph their reactor; we found it a much better 
model than their tame ducks, which refused to 
hold still for our cameras. (Their cows weren't 
very cooperative, either.) 
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5401B/9100 


5403A 


5402A 


5406B 


Multichannel-Analyzer System Capabilities 


MCA + 


MCA/Calculator 


MCA/ BASIC 


Multiparameter 




Calculator 


System 


System 


Analyzer 


UAIA ACCUMULATION 


| 









Number of analog-to-digltal conveners 


1 


1 


1 


1-16 


Maximum resolution (channels) 


B192 


8192 


8192 


24,576 


Data accumulation modes: 










Single parameter 


X 


X 


X 


x 


Multiparameter 








X 


Multiplex 








X 


Multichannel scaling 


X 


X 


X 


X 


Sampled voltage analysis 


X 


X 


X 


x 


DATA ANALYSIS 










Programming language 


calculator 


calculator 


BASIC 


EXECUTIVE 


Program storage media 
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Program console 


calculator 


calculator 


teleprinter 


teleprinter 


SpM d 


slowest 






- fastest 


Software expandability {subroutine language) 


calculator 


calculator 


BASIC, Assembly 


FORTRAN, ALGOL, 








Language 


Assembly Language 


SYSTEM CONTROL 










Analyzer functions programmable? 


No 


Yes 


Yes 


Yes 


DATA INPUT/OUTPUT 










Calculator peripherals 


x 


x 






Computer peripherals 






X 


x 


Analyzer peripherals 


x 


X 


x 




Coupler/controller peripherals 




X 






PRICE RANGE 


515 300 


$19,600 


S27.950 


535 000 up 



The smallest and least costly of the four systems pairs 
the multichannel analyzer with the HP 9100A/B pro- 
grammable calculator. With this system data from the 
analyzer can be entered directly into the calculator, either 
under control of a calculator program or at the user's 
command from the calculator keyboard. After the cal- 
culator has performed numerical analysis on the data 
to extract the useful information, the results can be dis- 
played on the calculator's cathode-ray tube or printed 
automatically by the calculator printer. A plotter can 
be added to the system to make permanent copies of the 
basic analyzer data. 

Next come two new and more powerful moderate- 
cost systems. In the HP 5403A MCA/Calculator Sys- 
tem the programmability and computing power are again 
provided by HP 9100-serics calculators, but instead of 
the one-way interface of the smaller system, there is a 
two-way interface through the HP 2570A or 2575A 
Coupler/ Controller. The interface gives the calculator 
random access to the analyzer's memory and control of 
important analyzer functions, and as a result the system 
is capable of a much higher degree of automation than 
the smaller calculator system. Also, the range of periph- 
erals and other instruments that can be added is much 
wider; besides the calculator peripherals most digital de- 
vices than can interface with the coupler/ controller can 
be made part of the system. For example, a teleprinter 
might be used to produce formatted reports. 

Both of the calculator systems are programmed from 



the calculator keyboard, which is relatively simple to 
use and requires no previous programming experience. 
Programs are stored on wallet-size magnetic cards for 
future use. 

In the other new system, the HP 5402A MCA/BASIC 
System, programmability and computing power are pro- 
vided by a general-purpose digital computer, the HP 
2114. The system is programmed in BASIC, a powerful 
but easy-to-learn conversational language. Users can 
easily program complex routines without becoming com- 
puter experts. Interfaces for all standard computer 
peripherals plug into the computer mainframe. The 
analyzer/computer interface is two-way. allowing the 
computer random access to the analyzer's memory and 
remote control of the analyzer. Thus the system is 
capable of a high degree of automation. 

The most powerful of the four types of systems is the 
HP 5406B Multiparameter Analyzer System. This 
computer-based system has a modular hardware and 
software architecture which makes it versatile and effi- 
cient. As a multiparameter analyzer it can include as 
many as four analog-to-digital converters like the one 
used in the 540 IB Multichannel Analyzer. It also has 
a multiplex mode in which as many as sixteen analog- 
to-digital converters can be interfaced with the computer. 
Its software, based upon a unique EXECUTIVE sys- 
tem, can easily be expanded by the user to include his 
own programs. 
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Identifying Radioactive 
Materials 

How multichannel analyzers work and why they 
need on-line data reduction 

Radioactive isotopes may emit many different kinds of 
radiation, but the kind most useful in nuclear studies is 
gamma-ray radiation A gamma-ray photon is emitted by 
the nucleus of a radioactive atom when a quantum transi- 
tion takes place between two energy levels of the nucleus. 
The energy of the photon is the difference between the 
energy levels involved in the transition. It is a fixed energy, 
an invariant of nature. A given radioisotope may emit such 
photons at one or more of these fixed energies, but no 
two isotopes emit the same combination of energies. Thus 
each isotope has a distinct gamma-ray energy spectrum, 
like a nuclear 'fingerprint,' which can be used to identify 
that isotope If It is present in a sample. 

Here's how a multichannel analyzer is used for identify- 
ing radioisotopes. The gamma rays emitted by a radioactive 
sample are sensed by a detector, The detector is a semi- 
conductor diode which emits an amount of charge pro- 
portional to the energy of the incident gamma-ray. These 
solid-state detectors are sometimes cooled to cryogenic 
temperatures to prevent changes in their structure and 
to reduce noise. The small amount of charge out of the 
detector is fed to a charge-sensitive preamplifier which pro- 
duces a voltage pulse whose height is proportional to the 
charge. After further amplification and shaping the pulses 
go to the multichannel analyzer, which measures the height 
of each pulse, converts it to digital form, and adds one 
count in whichever memory location, or 'channel', corre- 
sponds to that pulse height, or photon energy. 

The analyzer, in effect, grades' the pulses according to 
height in much the same way that a truckload of oranges 
might be graded according to size. If the oranges were 
allowed to roll down an incline over a series of bins, such 















1 

































































































F- "1 













t f t 



662 keV 1.17 1.33 
Cs 137 MeV MeV 



that the opening of each bin were a little larger than the 
opening of the preceding bin. then the smallest oranges 
would fall into the first bin, those a little larger into the 
next bin, and so on. One could then plot the number of 
oranges in each bin as a function of bin number or orange 
size. 

In analogous fashion, after counting pulses for a period 
of time ranging from minutes to days, depending upon the 
activity of the sample, the multichannel analyzer displays 
the gamma-ray spectrum of the sample, which is a plot of 
total counts versus channel number or pulse height or pho- 
ton energy. A typical spectrum has one or more peaks In 
it: the peaks are located at the characteristic energies of 
the photons emitted by the isotopes in the sample. If an 
ideal measurement were possible (it isn't), each peak would 
be very narrow — in fact, it would have no width at all. In 
real life, limitations of the detector and noise in the elec- 
tronics, mainly the preamplifier, tend to spread out the 
peaks into bell-shaped or Gaussian distributions. 

A measure of the analyzer's ability to separate closely- 
spaced peaks Is its digital resolution, or the number of 
discrete channels available to cover the range of energies 
expected in an experiment. The more channels an analyzer 
has, the higher is the analyzer's digital resolution. High- 
resolution analyzers like the HP 5401B have as many as 
8192 channels and produce results like the typical spec- 
trum shown here. The peak at 662 keV identifies cesium 
137. The two peaks at 1.17 and 1.33 MeV identify cobalt 
60. Although the Gaussian shape of the peaks isn't appar- 
ent here, the display could be expanded to show it. 

Once the analyzer has finished accumulating data, the 
mean value or centroid of each peak must be calculated 
to determine the energies of the photons emitted, and the 
total number of counts under each peak must be found to 
determine how much of each isotope is present. The latter 
calculation must be normalized to compensate for varying 
detector efficiencies and geometries. Thus two apparently 
simple calculations yield the peak locations and the rela- 
tive activity at each peak energy. 

A closer look at these calculations applied to an actual 
example shows that the situation is somewhat more com- 
plex. Background activity is always present, and this is 
accumulated along with the data from the sample. The 
background data may distort the shape of the peak and 
introduce error into the calculations. To obtain an accu- 
rate indication of the activity, the net area of the peak 
alone must be found. One technique for doing this assumes 
a straight-line function for the background and subtracts 
its area from the total, giving the corrected net area. A 
corrected centroid is calculated using a new peak found 
by subtracting the background data from the original peak 
data 

Normalization calculations can also become complex. 
Detector efficiencies change with energy, so an entire spec- 
trum may need amplitude normalization as a function of 
energy. 

Many more examples of data reduction techniques, such 
as curve-fitting, normalizing data for nonlinearities in the 
electronics, peak-to-background (signal-to-noise) enhance- 
ment, and so on. appear in current nuclear literature. Many 
of these calculations can be done quickly and efficiently 
by an on-line computer or calculator. This is the reason 
the systems described in this issue were developed. 
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The Systems In Action 

In the following paragraphs each of the four types of 
multichannel-analyzer systems is described in more de- 
tail, and examples are given to illustrate what each sys- 
tem can do. 

Analyzer/Calculator System 

In the low-cost analyzer/calculator system, a plug-in 
interface (HP 10667A) makes it possible to read ana- 
lyzer data into the calculator for processing. Data trans- 
fer can be initiated by FORMAT commands from the 
calculator keyboard or program. All calculator periph- 
erals can be used with the system. 

A typical application for an analyzer/calculator sys- 
tem might be found at a nuclear power station. Due in 
part to increased public concern, tighter standards are 
being placed on the radioactivity of the effluent from 
nuclear power plants. The amount of radioactive mate- 
rial that can be released into the environment from such 
plants is severely restricted, so much so that the most 
sensitive instruments available are needed to detect and 
measure any pollution and identify the isotopes respon- 
sible. Such instruments include high-resolution solid- 
state detectors and multichannel analyzers, combined 
with sophisticated data reduction equipment. 

To make certain that radiation limits are not exceeded 
routine monitoring of gaseous and liquid effluents is 
necessary. Samples are taken and analyzed for their 
radioactive constituents by measuring the various types 
of radiation being emitted. Approximately 50% of the 
radionuclides of primary concern emit gamma rays. Each 
sample analysis for these isotopes results in a gamma-ray 
spectrum on the multichannel analyzer. The analyzer 
data must then be reduced to determine how much of 
each radionuclide is present. The efficiency and geometry 
factors of the radiation detector and the data accumula- 
tion time arc known factors which arc used in the data 
reduction. The number of counts in a peak in the spec- 
trum, corresponding to the intensity of the gamma-ray 
radiation at a particular energy from a known isotope, 
must be converted to suitable units. This information, 
along with an experiment number, a sample identification 
number, the time, and detector constants, must be re- 
corded in a standard format. 

A simple calculator program recorded on a magnetic 
card can be used over and over again for both the data 
reduction and the data recording. Fig. I shows such a 
program. The user simply inserts the program into the 
calculator via the magnetic card, and enters the proper 



constants and identification numbers into the calculator 
via the keyboard. The result is a standard printout for 
each sample (Fig. 2) showing the required information. 

Another experimental record often desired is a nor- 
malized plot of the analyzer spectrum, such as the log 
or square root of the accumulated analyzer data. This 
plot can be recorded by the HP 9 125 A Calculator Plot- 
ter or its higher-speed relative, the 9I25B. 
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Fig. 1. Calculator program tor calculating isotope radio- 
activity trom the gamma-ray spectrum accumulated by a 
multichannel analyzer. This program is for the simplest 
possible analyzer/calculator system, which consists of 
analyzer, calculator, and interface kit. 
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Multichannel Analyzers and People 


Multichannel pulse-height analyzers were developed for 


■ 


Radiation-shield design 


nuclear research and their principal uses are still in the 


■ 


Reactor physics 


nuclear field. However, they affect our lives in many other 


■ 


Detector design and test 


ways, too. Here's a partial list of their applications. 


■ 


Nondispersive X-ray fluorescence 


■ Radionuclide assay 




Metallurgical surface physics 


Monitoring nuclear plant effluent 


■ 


Mossbauer spectroscopy 


Design, test, and construction of nuclear reactor 


■ 


Magnetic tape testing 


cores 


■ 


Noise and drift measurements in electronic circuits 


fuel-element leakage testing 


■ 


Pulsed-radar return signal analysis 


fuel-element scanning to find power-density dis- 


■ 


Particle size analysis (used in pollution studies) 


tribution in reactor core 


■ 


Medical sclntiphotography (organ imaging) 


neacior ruei reprocessing 






Nuclear materials management and inventory con- 




Nuclear physics 


trol 




isotope half-life determination 


Nuclear waste disposal control 




Nuclear medicine 


Environmental monitoring and fallout detection 




radioisotope dilution curves 


Medical diagnosis 




hepatic and renal blood-flow measurements 


whole body counting 




liver-function test (external hepatogram) 


radiopharmaceutical assay 




kidney-function test (renogram) 


Neutron activation analysis (used extensively in pol- 


■ 


Prospecting 


lution studies) 


■ 


Process control 


■ Nuclear physics research 


■ 


Communications-link analysis 


■ High-altitude physics research 


■ 


Optical transmission analysis. 



HP 5403A MCA/Calculator System 

The HP 5403 A System consists of an HP Series 9100 
Calculator, an HP 2575A or HP 2570A Coupler, an 
HP 540 IB Multichannel Analyzer, and interface boards 
which plug into the coupler and connect to the calcu- 
lator and the multichannel analyzer. In addition, periph- 
erals such as teleprinters, high speed paper tape punches 
and tape readers, digital printers or additional multi- 
channel analyzers can be added to the system by means 
of interface boards that plug into the coupler. Calculator 
peripherals such as the HP 9120A Printer, the HP 
9I25A/B Plotter and the HP 9 101 A Extended .Memory 
are compatible with the system, as are existing peripheral 
devices which connect directly to the multichannel ana- 
lyzer. Fig. 3 is a photograph of a 5403A system. 

All of the components in this system can be con- 
trolled from the calculator program or keyboard via 
FORMAT command sequences, or the PRINT com- 
mand in the case of the HP 9120A Printer. Multichan- 
nel analyzer functions which can be controlled from the 
calculator include the ERASE, ACCUMULATE. DIS- 
PLAY. Serial Output. Parallel Output. Serial Input, and 
TRANSFER modes, and START and STOP. In addi- 
tion, the calculator can access any memory channel in 
the analyzer and read or write data into the channel, 
can transfer data between the analyzer and a coupler 
peripheral or the calculator, and can check the operat- 

6 
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Several often-used nuclear data reduction programs 
are available lor the analyzer/calculator system. These 
include routines which calculate the net area under spec- 
tral peaks, the full width at half the maximum height 
(FWHM) of the peak, and the centroid of the peak. 
Many other routines can be written by the user, of 
course. For very long programs the HP 9 101 A Extended 
Memory may be needed. It increases the number of stor- 
age registers in the calculator system by 248; this is 
enough to store 248x14=3472 program steps. 




Fig. 2. Printout trom program of Fig. 1. 




Fig. 3. HP 5403A MCA/ Calculator System consists ot an- 
alyzer, coupler/ controller, and calculator. This system is 
capable ol a higher degree ol automation than the mini- 
mal analyzer I calculator system. The two-way coupler 
permits calculator control ol analyzer functions as well 
as transfers ot data among the analyzer, the calculator, 
and peripheral devices. 

ing status of the multichannel analyzer. 

As an example of what the system can do, take the 
nuclear-power-station monitoring problem again. With 
the small analyzer/calculator system the analyzer had 
to be set up manually by the user. The calculator pro- 
gram had to be rerun for each peak in the resulting spec- 
trum, and the user had to enter the time, the detector 
constants, and other parameters manually. 

The 5403A system, on the other hand, can do the 
entire job automatically, with the calculator program 
controlling the data accumulation, processing, and re- 
cording. Little training is required to operate the sys- 
tem once the program is written, since the program can 
be stored on wallet-size magnetic cards. 

In the typical situation, repeated measurements are 
required on the activity of effluent samples. The data 
accumulated must be analyzed to determine what the 
activity of several isotopes in the sample are. and the 
data must be recorded on a peripheral device in a suit- 
able form. If the activity of the source exceeds a certain 
value, an indicating device is to be triggered. 

Shown in Fig. 4 is a flow chart indicating how the 
5403A system would solve this problem. Once the in- 
struments arc set up and the program is entered into 
the calculator via a magnetic card the operation is auto- 
matic. The crystal-controlled timer in the analyzer can 
be preset from the calculator for any accumulation time 
up to 9.999 minutes, in 0.01 -minute increments. The 




Fig. 4. Flow chart ol calculator program lor analysis ol 
gamma-ray spectrum by HP 5403 A MCA/ Calculator Sys- 
tem. The calculator controls the analyzer, computes the 
activity ol all known isotopes in the sample, and triggers 
an alarm II the activity of any isotope is above a limit. 
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random-mcmory-access feature of the multichannel an- 
alyzer allows the calculator to take data from selected 
portions of the accumulated spectrum for data reduc- 
tion and analysis. The data is then recorded in a suitable 
format on the data recording device. If the data indicates 
isotope activity higher than is allowed, the calculator 
can command the coupler to trigger an alarm. 

There are many other applications for the 5403 A 
System, for on-line data reduction as well as automatic 
control of experiments. Several application programs 
have been written to help reduce analyzer data. One of 
these is for spectrum stripping. Here a known spectrum 
is multiplied by a constant for normalizing, and is sub- 
tracted point by point from a measured spectrum. The 
result is then written back into the analyzer memory in 
place of the measured spectrum. This type of data re- 
duction technique is useful in lower-resolution systems 
to identify the various isotopes represented in an un- 
known spectrum. Two address registers are provided 
on the analyzer/coupler interface for ease in program- 
ming such routines as spectrum stripping and integration. 

HP 5402A MCA/BASIC System 

The HP 5402A MCA/BASIC System combines the 
capabilities of the 540 IB Multichannel Analyzer with 
the versatility of a general -purpose digital computer. 
BASIC-languagc programming guarantees that the user 
does not have to become involved in how the computer 
operates. 

The HP 5402 A System (Fig. 5) consists of an HP 
540 IB Multichannel Analyzer, an HP 2114 Computer 
with 8192 words of core memory, a teleprinter, and 
interfaces. The computer and the analyzer can perform 
cither as stand-alone instruments, or together as a sys- 
tem with the computer controlling the analyzer and 
doing data reduction. For larger applications, several 
analyzers can be controlled by one computer. 

With the analyzer's FUNCTION switch in the EX- 
TERNAL position, major functions of the analyzer can 
be controlled by the computer. These include START 
and STOP, and the ERASE. ACCUMULATE. DIS- 
PLAY. Serial Output. Parallel Output. Serial Input, and 
TRANSFER modes. In addition, the computer can 
randomly access any channel of memory in the analyzer 
and transfer a word of data between that analyzer chan- 
nel and the computer, or transfer blocks of data between 
the analyzer and a computer peripheral. The operating 
status of the analyzer can also be monitored by the com- 
puter's BASIC program. 

HP Instrument BASIC is a computer language which 




Fig. 5. HP 5402A MCA/BASIC System consists of ana- 
lyzer, general-purpose digital computer, teleprinter, and 
interfaces. BASIC-language programming makes the sys- 
tem easy to operate. The analyzer and computer can be 
used separately or together as a system with the com- 
puter controlling the analyzer and doing data reduction. 

is well-suited as the man-machine interface in a multi- 
channel-analyzer system. It is simple to learn and use, 
yet powerful enough to perform nearly all desired com- 
putations. Analyzer functions are controlled by means 
of simple CALL program statements. 

The nuclear power examples discussed in connection 
with the calculator-based systems can also illustrate what 
the MCA/BASIC System can do. In addition to all the 
functions performed by the calculator in those examples, 
the computer-based system can perform more complex 
routines. For example, it can compare the calculated 
spectral peak energies with known values to find peaks 
that would otherwise be unrecognizable due to spectral 
noise or statistical fluctuations. Fig. 6 shows the BASIC 
program required to do this. The technique consists of 
correlating a characteristic spectral peak with the ex- 
perimental spectrum.' At points of high correlation a 
peak is said to exist. Fig. 7 shows a typical experimental 
spectrum and the correlation spectrum resulting from 
the BASIC program. Peaks in the correlation spectrum 
have the same mean values and proportional intensities 
(areas) as the original spectrum. The program computes 
these means and areas and prints them on the teleprinter. 

HP 5406B Multiparameter Analyzer System 

The HP 5406B is a modular computer-based single- 
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Fig. 6. B/4S/C program tor HP 5402A. The system accu- 
mulates a spectrum, then correlates it with the spectra of 
isotopes known to be present. This gets rid ot background 
data that might otherwise mask small peaks and make it 
impossible to measure their activities. 

parameter or multiparameter multichannel analyzer sys- 
tem. Its modular hardware and software give the user 
optional system configurations that include as many as 
sixteen analog-to-digital converters, a display subsystem, 
most computer peripherals including paper tape, mag- 
netic tape. disc, line printer, and others, and a variety 
of computer mainframe configurations with as many as 
24.576 words of data storage. Modular software makes 
it easy for the user to generate efficient real-time systems 
which may include programs supplied with the system, 
or programs written by the user, or both. Many often- 
used routines are supplied; examples are Gaussian curve 
fitting and multiparameter accumulate. User-written pro- 
grams can be in FORTRAN. ALGOL, or HP Assembly 
Language. 

Fig. 8 is a photograph of a 5406B system. All of its 
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Fig. 7. Analyzer displays ol portion ol raw spectrum (bot- 
tom) and same portion ol spectrum alter correlation by 
program ol Fig. 6. Correlation spectrum peaks have same 
mean values and same relative activities (areas) as peaks 
ol original spectrum. 

functions are controlled from the teleprinter keyboard 
by means of 5406 EXECUTIVE-Ianguage commands. 

An example of the use of a 5406B system is nuclear- 
reactor fuel-rod scanning. The rods have been used in 
an operating reactor core, and the analyzer system mea- 
sures their power densities as a function of position. The 
rods are extremely 'hot; so data rates are high. Four 
analog-to-digital converters are used; each views the rod 
from a different side. Using the analyzer's multiplex 
mode, the four 'views" are accumulated in four different 
segments of memory. Then the area of each peak in the 




Fig. 8. HP 5406B Multiparameter Analyzer is a powerful, 
modular, computer-based nuclear analysis system. It can 
include as many as 16 analog-to-digital converters and 
most computer peripherals. 
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resulting spectra is computed for each view as an indica- 
tion of power density. The results arc printed on the 
teleprinter and the original four spectra are stored on 
magnetic tape for future reference. 

Fig. 9 shows the 5406B EXECUTIVE program for 
this task. The operations MOVE and ANALYZE were 
written by the user in FORTRAN. Also shown in Fig. 
9 is the tabular format of the output. 



1- 



Mil 




Fig. 9. 5406B EXECUTIVE-language program and typical 
printout lor luel rod power density measurement using 
lour analog-to-digital converters. 
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Very High and Very Low Resistances - 
Why and How They are Measured 

These extremes of the resistance range can tell a great deal about properties 
of materials and the quality of semiconductor contacts. 

By Yoshihisa Kameoka and Jean E. Bonhomme 



The resistor, the simplest passive element, was first 
defined by an obscure German physicist. Georg Simon 
Ohm. His pamphlet, published in 1827, is entitled "Die 
galvanischc Kctte mathcmatisch bcarbeitet!' In it are 
the results of his efforts to measure currents and volt- 
ages, and to describe and relate them mathematically. 
One of his results is the fundamental relationship called 
Ohm's Law. 

The same result is said to have been discovered 46 
years earlier by a brilliant recluse. Henry Cavendish. 
But Cavendish's works were not discovered until long 
after both he and Ohm were dead. 

Ohm's law states simply that the voltage across var- 
ious types of conducting materials is directly propor- 
tional to the current flowing through the material. Stated 
another way. if the same potential difference is applied 
across a material, different currents will result. Thus, 
the relationship familiar to all electrical engineers: 

R=j-or (1) 

R = mL'i 'Q ' (2) 

Where R, the resistance, is the parameter responsible 
for the differences in current or voltage mentioned above. 
In Eq. 2, m — mass. L — length, t = time and Q = charge. 

Accurate and convenient resistance measuring instru- 
ments based upon Ohm's Law are in common use. Mea- 
suring resistance over a wide range, from less than one 
ohm to several megohms, is generally easy. Much more 
difficult have been applications where very low resist- 
ance, less than 1/100 ohm and very high resistance, a 
million megohms and higher, must be measured. Relay 
contact and semiconductor resistance are among the 
applications requiring low resistance measurements. In- 



sulating and dielectric materials, and high-voltage com- 
ponents are tested for high-resistance characteristics. 
Two Hewlett-Packard instruments, the HP Model 4328A 
Milliohmmeter and the HP Model 4329A High Resist- 
ance Meter aim at making measurements at these ex- 
treme ranges easier and more accurate than in the past. 

Low-Resistance Measurements 

The HP Model 4328A. Fig. 1. is a solid-state instru- 
ment with a range from 1 mf> to 100 CI. It uses the 
4-terminal (Kelvin) principle by which constant current 
is applied to the device under test, and the voltage drop 
is measured. A special probe is used, so that only two 
leads need be connected across the unknown. Errors 
due to contact potential and thermoelectric effects arc 
substantially eliminated by the use of a low level ac 
(1 kHz) measuring signal. Excessive voltage could 
change the surface characteristic of the sample under 
test. 

An IEC (International Electrotcchnical Commission) 
standard recommends that when contact resistance is 
measured by ac. the frequency should be 1 kHz ±200 
Hz. the rms value of current be 1 A or less, peak values 
of voltage be 20 mV or less, and instrument accuracy 
be ±10% or better. The Model 4328A meets these 
requirements. Voltage drop across the unknown is less 
than 200 ,«V peak on all ranges and does not exceed 
20 mV under any condition. Power dissipation in the 
devices under test is so low (23 i.W max) that the instru- 
ment is commonly used to check explosive devices such 
as fuses and squibs. 

A phase detector circuit is used in the Model 4328A 
which detects only the voltage drop in phase with the 
current through the device, and rejects the reactive com- 
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Fig. 1. This Hewlett-Packard Model 4328A Milliohmrneter 
measures low resistances from 1 milliohm to 100 ohms 
lull scale with a resolution ol 20 iiohms. 



ponent. Thus accurate measurements may be made of 
unknown resistances having reactance up to twice the 
full scale of the range in use. 

Principle of Operation 

The Model 4328A may be divided into two major 
sections. Fig. 2. A 1 kHz oscillator section serving as 
a constant current source to the sample, and a volt- 
meter section with a phase-discriminating circuit which 
amplifies the voltage drop across the sample, synchro- 
nously detects it. and deflects the ohm-calibrated meter. 

Assume that an unknown with an impedance Z t is 
connected between the current terminals of the probes. 
A closed loop is formed by the reference resistance R, 
of the instrument. Z, and transformer T2 secondary 
winding. A positive feedback loop is formed as both 
ends of R, are connected to the input terminal of high 
impedance amplifier, the power amplifier of the next 
stage and output transformer Tl. The circuit starts to 
oscillate at 1 kHz ± 100 Hz. the frequency determined 
by capacitance C„ and inductance L„ of Tl primary. 

The value of R. is selected by the range switch so as 
to be 500 times the value of the selected range. As a 
result the oscillator is a constant current source to Z s . 
Voltage across the reference resistance R, is always 
100 mV peak, thus when the value of Z, is equal to the 
full-scale value, the voltage applied to Z, is 200 mV peak. 
The input impedance of the high-input-impedance am- 
plifier is sufficiently high so that voltage changes across 
R, become negligible. If the probe leads are open, oscil- 
lation stops to prevent excessive voltage across the 
sample, and in the absence of oscillations, a Schmitt 
trigger in the circuit causes the meter to read towards 
infinity (full scale). 



The 20 mV peak limiter, connected to the tip of the 
current terminals, limits the ac voltage applied to the 
sample to 20 mV peak even if the range setting is incor- 
rect. This is done by comparing the 1 kHz output level 
to a reference diode. When a certain level is exceeded, 
the signal is amplified, rectified and the resulting dc level 
is used to reduce the oscillator gain. 

The voltage drop across Z, is applied to the ac am- 
plifier via transformer T2, and the amplified signal fed 
to one input of the phase detector. The other input is 
connected to a flip-flop driven by the 1 kHz oscillator. 
Because the output of the ac amplifier is synchronously 
rectified, the detected voltage is in phase with the cur- 
rent flowing through the reference resistance R,. 

An unknown sample can be defined as Z, = R I + jX t . 
and the resulting voltage drop V across the sample will 
be V=V K -\- jV r . But the phase detector detects only 
the voltage drop V,,- which is in phase with the current 
through R,. Therefore, errors due to jX r will be less 
than 1 % as long as its reactance value at 1 kHz is not 
more than twice the set range of the Model 43 28 A. 

Some Applications 

A typical circuit for measuring dynamic resistance 
of various semiconductor devices is shown in Fig. 3. 
Devices commonly measured include switching diodes, 
power diodes and point contact diodes. Resistivity of a 
semiconductor wafer is easily measured using a four- 
point probe. 

High-Resistance Measurements 

Related to resistance is resistivity p, which is the 
characteristic of a material rather than of a particular 
specimen of that material. The resistivity of copper, for 
example, is 1.7 X 10'' ohm-cm. and that of quartz is 
about 10" ohm-cm. Resistivity is of importance when 
considering the fundamental behavior of materials, as in 
solid-state physics. It is also important when investigat- 
ing the interior behavior of irregularly-shaped objects. 

Surface and volume resistivities of solids such as 
printed circuit board materials, plastics, resins, refrac- 
tories, semiconductor materials, oils, solvents and fabrics 
require high-resistance measurement capability. Other 
more common applications include measurements of 
capacitor insulation, and relay coil resistance. Accurate 
high-resistance measurements arc needed to examine 
effects of temperature, humidity, and surface conditions 
on a variety of insulating devices and materials. More 
important, most high-resistance materials have non- 
linear characteristics, i.e. voltage versus resistance. Thus 
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Fig. 2. Functional block diagram 
of the Model 4328A shows how 
a feedback loop is formed when 
Z t is plugged in. Inductance of 
the primary winding of T1 and 
C 0 determine the frequency of 
oscillation. The modulator pre- 
vents the meter voltage (V r ) 
from varying in the presence of 
small changes in teedback loop 
or oscillator output. 




Fig. 3. A typical circuit used to measure the dynamic re- 
sistance of a diode. 

voltage measurement is very important. 

Measuring of high resistance is essentially a problem 
in measuring extremely small currents. Until recently 
the vacuum-tube electrometer, the vibrating reed capaci- 
tor and the field effect transistor have been the only 
devices available to provide the high input impedances 
required. 

Further, it has been necessary to design dc to ac con- 
verters to provide ac amplification of the outputs of 
tubes and transistors used in electrometers. These con- 
verters (usually choppers) are required to eliminate the 
errors inherent in dc amplification of such low-level 
signals. The combined inherent instabilities, drift, leak- 
age currents, grid currents and converter problems have 
justified devising an entirely new principle in this type 
of measurement. 

As used in the HP 4329A High-Resistance Meter. 
Fig. 4. this principle is called the 'self-oscillating para- 



metric impedance converter! Fig. 5. This circuit provides 
both a high-input impedance and an ac output available 
for stable, drift free amplification. A variable-capaci- 
tance diode (varicap) bridge serves as the parametric 
element. 

These diodes are variable capacitance diodes designed 
to have a known junction capacity versus voltage char- 
acteristic variable with applied voltage. Normally this 
characteristic is non-linear and increases with increas- 
ing forward bias. Here the diodes are used as variable- 
capacitors since the applied forward bias does not extend 
into the conduction region of the diode. 

Fig. 6 illustrates the diode characteristic curve show- 
ing the capacity of the depletion layer or 'operating 
region' as a function of applied diode voltage. 

Principle of the Parametric Converter 

A parametric amplifier is defined as a network capa- 
ble of transferring energy from a high frequency source 
to an applied signal at a lower frequency by virtue of 
a periodically varying energy storage element. The time 
varying element need not be linear. 

The Model 43 29 A parametric converter circuit can 
be simplified to a basic bridge circuit. Fig. 7. CI and 
C2 are diode capacitors; the two resistance arms are 
equal. The bridge unbalance voltage is the input to an 
amplifier, the amplifiers output is the 'supply' for the 
bridge. 

To understand the parametric converter, consider the 
result of making a very small change in the value of C T . 
If for example. C T is raised in value, the net C in the 
lower arm now is greater than that of CI. Thus the 
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Fig. 4. High resistances to 2x10 ' ohms are measured 
with the Hewlett-Packard 4329A High Resistance Meter. 
Automatic range and multiplier indicators are lighted to 
make them easy to read. 




Fig. 5. Basic principle ol the measuring technique used 
lor high resistance measurements. 



voltage across the lower arm will tend to decrease to a 
value less than that across CI. Point 1 tends to become 
more negative. This represents a bridge unbalance, and 
the amplifier sees this as a net negative-going signal at 
its input. 

This signal is inverted and amplified and the bridge 
"supply voltage' changes. This change in the net voltage 
between point 3 and 4 is seen by CI as an increase in 
forward bias, and by C2 as a decrease in forward bias. 
CI therefore tends to increase its capacitance, and C2 
tends to decrease its capacitance. Thus the effect of the 
fed-back voltage is to cancel the initial unbalance due 
to the change in C T . 

However, the circuit constants (R. L, C) determine 
the speed at which the response can catch up with the 
initial unbalance, and by the time the bridge is re- 




Fig. 6. Characteristic ol the variable capacitance diode 
used in parametric converter. The diode used in this 
instrument has a Co ol 100 pF. 



balanced and quiescent, the total C value of the upper 
arm is greater than that of the lower and the cycle re- 
peats with all polarities reversed. This action continues 
as a normal oscillator, at a frequency and amplitude 
determined by the initial change in C T and by other 
component values. This type of oscillation is called a 
parametric oscillation, in that the modulation of the 
reactance parameter (C) determines the behavior of the 
oscillator. 

In applying this concept to a practical high-resistance 
measuring circuit, only a few simple changes are neces- 
sary. Point 1. the junction between C' an d C2 can be 
considered the input point, where a voltage to be mea- 
sured is applied to the circuit. 

Any change in voltage at this point will be seen as a 
further unbalance in the bridge. This additional unbal- 
ance will increase the signal to the amplifier and result 




Fig. 7. Basic parametric converter bridge circuit. 
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in an increase in the amplitude of oscillations. 

It remains only to develop a dc voltage proportional 
to the increase in oscillation amplitude, suitable to drive 
a meter circuit. The high-impedance input is developed 
by the non-conducting variable-capacitance diode, aided 
in a practical circuit by suitable isolating capacitors, high 
resistance circuit mountings, and augmented by nega- 
tive feedback. The stable gain characteristic of the Model 
4329A is accomplished in the standard ac feedback- 
controlled amplifier manner, made possible by convert- 
ing the dc input current to a change in the amplitude of 
converter oscillations. The actual circuit used in the HP 
4329A. Fig. 8. operates like the analogous circuit. 

Practical Circuit 

The converter consists of a transformer bridge. Diodes 
CR1 and CR2 are two arms, and the two halves of 
center-tapped Tl primary winding are remaining two 
arms. C6 and C7 block dc from the transformer wind- 
ings. C T is the variable capacitor providing a slight un- 
balance to the bridge and is used to set the 'no-input" 
oscillation amplitude. Bridge unbalance is taken from 
the junction of C T and C„ to the oscillator/amplifier cir- 
cuit through series resonant circuit C„/L,.. which sets 
the oscillation frequency at approximately 100 kHz. 

Output of the amplifier is fed back to the bridge by 
transformer action (from M3 to M2. Ml). Oscillations 
are also transferred by transformer action to the input 
of meter amplifier circuit rectified, and the resultant 
direct current drives the meter. The meter circuit cur- 
rent is set such that 200 mV at the converter input gives 
full-scale deflection. 



DC negative feedback through the RC network to 
the pair of 2 Mn resistors stabilizes the circuit against 
variations in circuit parameters, particularly the variable- 
capacitance diode characteristics. DC feedback to the 
guard at the input point adds to the total effective insu- 
lation resistance of the insulator which isolates the con- 
verter input lead from the leakage resistance of the 
printed circuit board. 

This feedback equals the input voltage itself, and is 
applied to the conducting surface around the periphery 
of the insulating cylinder containing the input junction 
point. Since the voltage at the inner surface of the cyl- 
inder is equal to that at the outer surface, no leakage 
occurs through the teflon cylinder and its insulating 
qualities are enhanced. 

Covering a resistance range of 500 ki> to 2 .:< 10" ; <}, 
the Model 4329A. also has seven discrete test voltages 
from 10 to 1000 V. Test voltages are selected by a 
front-panel switch. 

Resistivity Cell 

Volume and surface resistivities can be measured by 
the Model 4329A using a standard mercury cell. An ac- 
cessory Model 1 6008 A Resistivity Cell is more conven- 
ient, although a little less accurate. The cell's main elec- 
trode is coated with conducting plastic and is surrounded 
by a guard ring. Spring loading the upper electrode pro- 
vides pressure to measure unknown samples up to 7 mm 
thick, and a switch selects volume or surface measure- 
ment. Contact between the conductive plastic and the 
surface of the device under test may not be as good as 
the mercury cell contact, however. 
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PARTIAL SPECIFICATIONS 

Model 4329A 
High Resistance Meter 

RESISTANCE RANGE: 500 k(J to 2 X 10" 0. 
CURRENT RANGE: 0.05 pA-20 *A. 

RECORDER OUTPUT: 0 10 100 mV dc. proportional lo meter deflection. 

output Impedance 1 kO. 
POWER: 115/230 V, =10%. 50 or 60 Hz. 3 W 
DIMENSIONS: 6.5 In high X 7.8 In wide X 8 In deep 

(186 X 19B X 204 mm). 
WEIGHT: 7.7 lbs |3.5 kg). 

PRICE: Model 4329A High Resistance Meter $750.00 
ACCESSORIES FURNISHED 

16117A Low Noise Test Leads, composed of Insulated BNC to Alli- 
gator clip lead, and banana plug to Alligator clip lead. On BNC 
lead, outer conductor comprises guard shield. Guard isolated from 
low resistance alligator clip sleeve. 

Model 16008A 
Resistivity Cell 

RANGE: For sample thicknesses approx 1 mm, up lo 4 X 10 * li — cm In 
volume resistivity. Up to 4 X 10" t! In surface resistivity. 

TEST LEADS: Provided; cell top cover Inlertocks with test voltage. 

MAIN ELECTRODE: 50 mm diameter, conductive-plastic coated. 

UPPER ELECTRODE: 100 X 120 mm. 

GUARD RING: 70 mm diameter. 

MAX TEST VOLTAGE: 1000 V. 
DIMENSIONS: 1.9 in high X 7.7 In wide X 6 In deep 
(48 X 195 x 152 mm). 

WEIGHT: 3.7 lbs (1,8 kg). 

PRICE: $200.00 

MANUFACTURING DIVISION: YOKOGAWA-HEWLETT-PACKARD, LTD. 
9-1, Takakura-cho 
Hachloll-shl 
Tokyo 192, Japan 

PARTIAL SPECIFICATIONS 

Model 4328A 
Milllohmmeter 

RANGE: 0.001 to 100 ohms full scale In a 1, 3, 10 sequence. 

RESOLUTION: Max resolution 0 02 on 0-1 scale, 0.05 on 0-3 scale. 

ACCURACY: ^ 2% of full scale. Accuracy is unaffected by series react- 
ance of samples with magnitude of up to 2 times lull-scale resistance 
value. 

MEASURING FREQUENCY: 1000 Hi =100 Hz. 
VOLTAGE ACROSS SAMPLE: 200 ,.V peak at full scale. 
MAXIMUM VOLTAGE ACROSS SAMPLE: 20 mV peak in any case. 

SUPERIMPOSED DC: 150 V dc max may be superimposed on samples 
from an external source. 

RECORDER OUTPUT: 0.1 V dc output at full scale meter deflection. 
Output resistance approx. 1 kO. 

POWER REQUIREMENT: 115/230 V =10%, 50 to 60 Hz. 1.5 W. 

PRICE: Model 4328A S495.00. Option 001, rechargeable battery opera- 
tion, add $25.00. 

MANUFACTURING DIVISION: YOKOGAWA-HEWLETT-PACKARD, LTD. 
9-1, Takakura-cho 
Hachloll-shl 
Tokyo 192, Japan 
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